Singular Points of Ordinary Differential Equations 

Lecture 7 Regular Singular Points of Ordinary 
Differential Equations 

Let us consider the Bessel equation 

~^>- < 6 - 22 > 

Let us seek the solution of the Bessel equation which is some kind of expansion around the point 
x = 0. 

When x 2 is much smaller than p 2 , we make the approximation of neglecting x 2 y as it is much 
smaller than p 2 y. In this approximation we set the right-side of (6.22) to zero. We note that the 
resulting ODE is equidimensional. Thus the solution of the resulting equation is of the form x s . 
Substituting this solution into the differential equation, we get 

s(s- 1) +s-p 2 = 0, 

or 

s = p, or -p. (6.23) 

Thus the approximate solutions are 

y(x) « x p oxx-p, (6.24) 

which are good approximations when x is small. 

Let us seek an improvement of the approximate solutions of (6.24). 
We express y as 

y(x) =J] a n x n+s . (6.26) 



As before, we shall understand that 

a_i = a_2 = ■■■ =0. 

We shall require that 

ao * 0. 

This means that aox s is by assumption the first non-vanishing term in the series solution of the 
Bessel equation. The series in (6.26) is called the Frobenius series. We shall show that, as we 
substitute (6.26) into the Bessel equation, we determine the index s as well as the coefficients a n , 
n = 1,2 • • •. The values of the coefficients a n ,n = 1,2 • • •, depend on the value of s as well as ao. 
There are two possible values of s, and hence we obtain, except in some special cases, two Frobenius 
series solution of the Bessel equation. The general solution of the Bessel equation is a linear 
superposition of these two solutions. This method is known as that of Frobenius. 

Problem for the Reader: 
Find the recurrence formula for the Bessel equation. 

Answer 

Determining the dimensions of the terms in the equation will help to reduce the chore of 
calculation. This is because as we apply two operators of the same dimension to a Frobenius series, 
the two resulting series can be added together without re-assigning the index of summation, as both 
operators change the power of x n+s by the same amount. Thus it will save both space and labor if we 
handle the operators with the same dimension at the same time. 

The only term in the Bessel operator (x 2 D 2 + xD + x 2 - p 2 ) of dimension 2 is x 2 . We have 

=2 a n X " +S+2 =2 a n-2X" +S - 
n n 

The rest of the terms in the Bessel operator are of dimension zero, and we will take care of all of 
them at the same step. We have 

{x 2 D 2 + xD -p 2 )y = ^[(n + s)(n + s-l) + (n + s) - p 2 ]a n x n+s ~ 2 . 

Thus the recurrence formula for the Bessel equation is 

(n + s + p)(n + s - p)a n = -a n -2- (6.27) 
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Singular Points of Ordinary Differential Equations 
We note that there are only two coefficients in (6.27): a n and a n -i ■ This is because there are only 
two different dimensions for the operators in the Bessel equation. 

To obtain the solution of the Bessel equation, we first set n = in (6.27). Since a_2 is zero and 
since ao * 0, we get 

(s+p)(s-p) = 0. (6.28) 
Equation (6.28) is called the indicial equation. It determines the index s to be either si or S2- We get 

si = p, s 2 = -p, 

in agreement with (6.23). Without loss of generality, we may assume that Rep > 0, and hence 
Resi > Res2- 

Next we write (6.27) as 

a n = %^ r,n > 0. (6.29) 

(n + s + p)(n + s -p) 

Problem for the Reader: 
Solve the recurrence formula (6.29). 

Answer 
Let us set n = 2m in (6.29). We get 

Q2(m-1) 

dim — .. s+ p s -p . ■ 

4(m + —j-)vn + -£-) 

By applying the above recurrence formula m times, we get with the help of (6.10) that 

_ ( yn a r( 1 + s/2 + p/2)F( 1 + s/2 - p/2) ( , - m 

2m 1 ; 4 m T(m + s/2 +p/2+ l)T(m + s/2 - p/2 + 1) ' 
The right-side of (6.30) depends on s, with the two values of s determined from the indicial 
equation (6.28) given by (6.23). 
Problem for the reader: 
Find a solution of the Bessel equation. 

Answer 
Setting s = p in (6.30), we get 

4T(m + p + l)r(m + 1) 
Thus one of the solutions of the Bessel equation is 

2m+p 



(f) 



M*) =S (-l) m r . /^, rr -pr , (6.32) 
r(m +/? + l)r(m + 1) 

where the immaterial constant multiple aor(l + p) has been ignored. The function J p (x) is known to 
be the Bessel function of order p. 

We may get a second solution of the Bessel equation by choosing s in (6.30) to be the other 
index —p. But there is a slightly easier way. The Bessel equation is unchanged as we replace p by -p. 
Thus another solution of the Bessel equation is obtained by replacing p in (6.32) by -p. The solution 
obtained is J- P (x). 

By (6.29), the ratio a n /a n -2 vanishes like lln 2 as n ->■ qo. Thus the series (6.32) converges for all 
values of x. Therefore, the series is not only useful as an approximate series when x is small. It is 
meaningful even if x is not small. However, when x is large, it is inefficient to calculate the values of 
the Bessel function J p (x) with the series (6.32). Also, the series hardly tells us anything about the 
qualitative behavior of the Bessel function for large values of x. To get a hold on the behavior of the 
Bessel function for large values of x, it is a good idea to obtain its expansion in the variable x~\ 
rather than the variable x. 

Consider now the general second-order linear homogeneous differential equation 

y" + c(x)y' + d(x)y = 0. (6.38) 

If c(x) and d(x) are both analytic at xo, xo is called an ordinary point of (6.38). And if either c(x) or 
d(x) or both have a singularity at xo, xo is called a singular point of (6.38). 

If xo is an ordinary point of (6.38), we may obtain the two independent solutions of (6.38) by 
seeking them in the form of a Taylor series expanded around xq: 
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y( x ) = ^ja„(x-xo)",a-\ = a -2 = ■■■ = 0. 

The coefficients a n ,n = 2, 3 • • •, are determined from the recurrence formula obtained by plugging 
the series into the differential equation. These coefficients are dependent on ao and a\, which are 
arbitrary constants. We get two independent solutions, the first being obtained with the choice of 
ao = 1 and a\ = 0, and the second one being obtained with the choice of ao = and a\ = 1. 

Let xo be a singular point of (6.38). There are two kinds of singular points. If xo is a singular 
point of (6.38) with (x - xo)c(x) and (x - xo) 2 d(x) both analytic at xo, then xo is called a regular 
singular point of (6.38). 

The Bessel equation (6.22) takes the form 

/ + ly + (i-£) y = o. 

Comparing with (6.38), we find that the coefficient c(x) and d(x) for the Bessel equation are x -1 and 
(1 - p 2 /x 2 ), respectively. The point x = is a regular singular point of the Bessel equation, as c(x) 
has a simple pole and d(x) has a double pole at this point. 

Retuning to (6.38). We may obtain the solutions of (6.38) expanded around a regular singular 
point xo in the same way we obtained the solutions of the Bessel equation expanded around x = 0. 
More precisely, we seek these solutions in the form of the Frobenius series 

y(x)=J^a n (x-x r +s , (6.39) 

n 

with 

a_i = a_2 = ■■■ =0. 
The series is convergent at least up to the next singular point of (6.38). 

Finally, if either (x - xo)c(x) or (x - xo) 2 d(x) or both are not analytic at xo, then xo is called an 
irregular singular point of (6.38). 
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